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ABSTRACT 

We study the effect of reddening on microlensed stars towards the LMC. If lenses 
are in the LMC disk, then the source stars should be at the far side or behind the LMC 
disk. Thus they should experience more reddening and extinction by dust in the LMC 
disk than typical stars in the immediate neighbouring lines of sight. We simulate this 
effect in a variety of models for the LMC stars and dust. We stress that the optical 
depth is not one number, but a function of the reddening of the survey stars. We 
discuss how these effects could be used to constrain the fraction of machos in the dark 
halo. The effect of patchiness of dust can be controled by working with faint stars in 
the smallest patches of the sky around the microlensed stars. This can be done most 
effectively with the Hubble Space Telescope in the ultra-violet. The non-detection of 
the reddening bias would be strongly in favor of machos in the Galactic halo. 

Subject headings: dust, extinction - Magellanic Clouds — Galaxy: structure 

1. Introduction 

One of the main puzzles of Galactic microlensing surveys is the poorly determined location 
of the lens population of the events towards the Magellanic Clouds. Currently there are two 
popular views on the issue: (a) the lenses are located in the halo, hence are likely baryonic dark 
matter candidates (Alcock et al. 1997); (b) both the lenses and sources are part of the Magellanic 
Clouds, hence are stars orbiting in the potential well of the Clouds (Sahu 1994, Wu 1994, Zhao 
1998a, b, 1999a, Weinberg 1999). These lead to drasticly different predictions the relative amount 
of baryonic dark matter in galaxy halos as compared to baryons in stars and gas of various phases 
in galaxies and intergalactic medium. Gould (1995) showed that the low dispersion of many LMC 
tracers limits the amount of self-lensing in a virialized disk of the LMC, as assumed in original 
models of Sahu and Wu. Zhao (1998a) suggested that this constraint can be circumvented by 
invoking a plausible amount of unviralized material, particularly tidally excited material in the 
immediate surrounding of the LMC disk, as either foreground lenses or background sources. Zhao 
(1998b) suggested that the polar stream seen in Kunkel et al.'s (1997) kinematic survey of the 
LMC carbon stars could contribute to lensing. Weinberg (1999) showed a specific model for 
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thickening the LMC disk by interaction with the Galaxy. Presently the observational status is 
still very confusing. Several observations suggest that our line of sight to the LMC passes through 
a 3-dimensional stellar distribution more extended in the line of sight than a simple thin disc of 
the LMC (Kunkel et al. 1997, Zaritsky et al. 1997, 1999 and references therein). Zaritsky et al. 
(1999) also include a review of the arguments and counter-arguments for these structures. 

To break the degeneracy of the models we need new observations which are sensitive to the 
thickness of the LMC and the location of the lenses. Several ways have been proposed, including 
studying the spatial distribution and magnitude distribution of the microlensed sources (Zhao 
1999a, Zhao, Graff & Guhathakurta 1999), the radial velocity distribution of the sources (Zhao 
1999b) and the reddening distribution of the sources (Zhao 1999c). Essentially if lensing is due 
to Galactic halo machos or a foreground object, the sample of lensed stars should be a random 
subsample of the observed stars in the LMC. On the other hand, if there are some tidal material 
behind the LMC disk, then the stars in the tidal debris can be strongly lensed by ordinary stars 
in the intervening LMC disk, hence the spatial and kinematical distributions of these source stars 
will be those of the tidal debris, rather than those of the LMC disk. It is suggested that the source 
stars in the tidal debris should stand out as outliers in the radial velocity distribution of the LMC 
(Zhao 1999b), or fainter-than-normal red clump stars in the color magnitude diagram (Zhao, Graff 
& Guhathakurta 1999). 

The idea about the reddening distribution is to measure the distribution of the reddening of 
individual LMC stars in small patchs of sky centered on the microlensed stars. Basicly some kind 
of "reddening parallaxes" can be derived for these stars from the line of sight depth effect, i.e., the 
dust layer in the LMC makes stars behind the layer systematicly redder than those in front of the 
layer. The method involves obtaining multi-band photometry and/or spectroscopy of fairly faint 
(19 — 21mag) stars during or well after microlensing. From these data we can infer the reddening 
of microlensed stars and neighbouring unlensed stars. The method, in some sense, is a variation 
of the distance effect suggested by Stanek (1995) for the Galactic bulge microlensing events. 

Here we extend the analysis of Zhao (1999c) and present simulations of the reddening of 
microlensed sources and show how reddening might be used to constrain the macho fraction in the 
halo. The models are described in §2, results are shown in §3. We summarize in §4 and discuss a 
few practical issues. 



2. Models 
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2.1. Density models for the dust, the stars and the lenses 

For the time being we shall assume a smoothed dust layer of the LMC with a sech 2 profile 
and a FWHM of w. Let Pd{D) be the dust density at distance D, then 

Pd (D)dD = Pd ( J D LM c)sech 2 ( D ~ dD, (1) 

where pd(DhMc) is the dust density at the mid-plane of the LMC. 

For the LMC stars, we assume a two-component model. Let plmc(-D) be the density of the 
LMC stars, then 

PLMC(£*) 

where 

L = 



Our LMC model includes a (thin) sech 2 -disk component of a surface density Sdisk an d FWHM 
thickness W, and an extra component with a surface density S C xtra and FWHM thickness Wgxtra- 
The extra component consists of two half-Gaussians, which are joined together at their peaks at 
D = -Dcxtra and truncated at the lower and upper ends at D = D m \ n and D = D max . This is 
intended to model any non-virialized material in the vicinity of the LMC (Zhao 1998a, Weinberg 
1999). To be most general, an offset between the extra component and the LMC disk, namely 
-Cextra / -^LMCj is allowed in our models. 

Let v*{D)dD be the number of LMC star in a line of sight distance bin (D, D + dD), then 
the LMC star number density will be given by 

v*(D)dD = C oPLM c(D)D 2 dD, (5) 

where Cq is a normalization constant. 

For the lens population, we include the contribution from the LMC stars and intervening 
machos in an isothermal halo. The density of lenses at distance D\ is 

Plens(A) = PLMC(A) + /machoPhalo(A), (6) 

where /macho is the fraction of machos in the halo, and Phaio(-D) is the dark halo density given by 

-l 

(7) 

The halo density corresponds to a surface density of dark matter (machos plus wimps) halo 
Shalo ~ 100M Q pc~ 2 towards the LMC. 
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2.2. Reddening vs. line of sight distance of a star and their distributions 

In general the reddening of a star is correlated with the distance to the star, with stars at 
the backside of the LMC seeing more dust and experiencing more reddening. But since the dust 
distribution is clumpy and the reddening distribution is patchy, any relation between the reddening 
and the line of sight distance is not a well-defined one-to-one relation, but has a significant amount 
of scatter due to patchiness of dust and measurement error. Nevertheless we can argue that stars 
in a small (e.g., 4" x 4") patch of sky will see the same set of dust clouds, so if we sort these 
neighbouring stars according their reddening, we also get a sorted list in distance. 



Stars in the LMC can be ranked according to an observable reddening indicator k ^ s . For any 
LMC star define 



where E(B — ^)lmc * s the observed reddening in B — V color towards the LMC star, after 
discounting the reddening by the Galactic foreground, and \E(B — 7)^/ is the average 



/LMC / 

over random stars in the small patch of sky centered on that LMC star, and is approximately 
the reddening of those stars at the mid-plane of the LMC. Surely the typical reddening 
(^E(B — ^)lmc) var i es widely from patches to patches on scales of arcmins (Harris et al. 1997 
and references therein), but the rescaled reddening K Q b s allows us to sort the stars not only in the 
same patch but also stars in different patches. A higher rank, i.e. a larger K bsj means a deeper 
penetration into the dust layer. Stars at the midplane have a rank k ^, s = 1, and those at the front 
or back side of the LMC disk have a rank K Q b s = or K b s = 2. 

Consider the reddening rank of stars in a smooth dust model. Let A^mc(D) be the dust 
absorption towards a LMC star at distance D (after discounting the absorption by the Galactic 
foreground dust), and let k(D) be the predicted reddening rank of the star, then 

K ( D ) = fWg) = Io%(D)dD (g) 
4mc(-Dlmc) J° LMC Pd (D)dD 

(cf eq. [l]), where Pd(DhMc) and ^lmc(-Dlmc) are the dust density and absorption at the mid-plane 
of the LMC. It is easy to show that 

k{D) = l + tanhf D ~f LMC ") (10) 



0.56w / 

= 0, D-D LUC <£w (11) 

= 1, D = D hMC (12) 

= 2, D-D LMC ^>w. (13) 

So the reddening rank k{D) increases from to 2 abruptly in the distance range -Dlmc ± f as a 
star enters and exits the dust layer. It turns out that this property is not limited to sech 2 profile 
dust, but generic for smooth dust layer of a symmetric profile. 
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Now to take into account of the scatter due to measurement error and any residual patchiness, 
we assume that the observed reddening has a simple Gaussian distribution centered on the 
predicted value from a smooth dust disk model. Then for a star at any distance D, the reddening 
rank k \, s is drawn from the following distribution function 

\2' 



B(n ohs ,D) 



1 
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exp 



(K ohs - k{ D)Y 

2a 2 



(14) 



which is a Gaussian with a constant dispersion a, and a mean n(D) as predicted from a smooth 
model. 

Now we integrate over LMC stars at all distance and bin them according to their reddening 
ftobs- Let -/V*(re b s ) be the relative frequencies of finding an unlensed star with a reddening K bs, 
then 

A^obs) = / dDv*(D)B( Kohs ,D). (15) 



2.3. Optical depth as a function of distance and reddening 

The microlensing optical depth t(D s ) of a source star at distance D s is defined by 

s Sty (16> 

where v*{D s )dD s and v s (D s )dD s are the numbers of stars and microlensed sources in the distance 
bin (D s , D s + dD s ) respectively (cf. eq. ||). It is well-known that in the standard full-macho 
isothermal halo model / mac ho = 1> Piens = /^haio (cf. eq. |6|), and the microlensing optical depth 
r std = 5 x 10 -7 (cf. e.g., Paczyhski 1986). It is convenient to rescale the optical depth of a model 
with this standard value. The microlensing optical depth t(D s ) is then given by 

t(d)=t J Ds dD lPlens (Di)(Ds-Di)Di/Ds (17) 

Std f MC ^ftalo(A)(OLMC-A)A/5LMc' 

The optical depth t(D s ) increases with the source distance. Suppose all lenses are at distance 
Di, then 

t(D s ) ~ {D s - Di). (18) 

Unfortunately, it is difficult to observe this effect because it is hard to measure the distance 
to a LMC star at the accuracy of 1 kpc. But fortunately the distance to a star is correlated with 
the reddening, and the reddening of a star can be measured accurately to 20%. So we have an 
interesting effect that the rate of microlensing is correlated with the reddening of the stars. For 
each bin in terms of the rank of reddening (K b s , K Q bs + d^obs) 5 we can define an optical depth 

/ . _ N s (K ohs ) _ J dD s u s (D s )B(K ohs ,D s ) 
r obB {n oh8 ) _ - $dD<D)B{KQ ^ D) , U»J 



- 6 - 



where iY*(fc b s ) and N s (K Q b s ) are the relative frequencies of finding an unlensed star and a 



microlensed source with reddening k q ^ s respectively (cf. eq. 15). 



Note r bs(^obs) is a function measurable observationally because we can always bin the 
microlensed events according their observable reddening, and get their frequencies iV*(K bs) and 
^Ys(ftobs)- Of particular interest is T o b s (0 < K Q bs < 1)- Note stars at the mid-plane of the LMC 
disk have average reddening, so K Q bs ~ 1 ± 3<r at 3c" confidence level. In comparison stars well in 
front of the dust layer have K Q bs < 3<r, and stars well behind the dust layer have K Q bs > 2 — 3a. 
So T o b s (0 < K Q bs < 1) is the optical depth of a thin layer of stars co-spatial with the dust layer 
and slightly closer to us the mid-plane of the LMC with -Dlmc — w/2 < D < -Dlmo Since the 
thickness of the dust layer w ~ 100 — 200pc is very small, these stars are virtually at the same 
distance. So we have 

7"obs(0 < Kobs < 1) ~ t(D LMC ) > /macho^std- (20) 

So we have turned the theoretical quantity t(-Dlmc) to an observable, and this observable optical 
depth r o bs(0 < k q ^ s < 1) sets an upper limit on the amount of machos in the foreground of the 
LMC. The above inequality reduces an equality in the absence of any LMC lenses in front of the 
LMC disk. 

If the picture that the LMC is a thin disk without any extra material in the front or back is 
correct, then we should find that r bs(K bs) is at a constant level « t(Dlmc) = /macho^std for all 
reddening. 

Presently the survey teams have not studied the reddening dependence of the optical depth, 
so they quote only the observed optical depth averaged over all LMC stars, which is only one 
number. Let this be (r b s ), then current data require 



(Tobs) ~ iS:.y D ^(D) ~ idK 0hs N*(K 0hs ) ~ u - 5Tstd ' 



where the factor 0.5 comes from the fact that the observed microlensing optical depth account for 
about one-half of that of the standard macho dark halo model r st d (Alcock et al. 1997). 

We would like to know how much of the observed optical depth is due to stellar lenses in the 
LMC vicinity and how much is due to machos in Galactic halo. If we define /* as the fraction of 
the observed microlensing depth due to stellar lenses, then the fraction of machos in the dark halo 
is given by 

/macho = (T ° bs) ^ ~ U) = 0.5(1 - /*). (22) 
Tstd 

So for a given set of stellar distribution parameters S ex tra etc, we can predict the optical depth 
due to stellar lenses /*r st( j, we can then determine /macho by fixing the overall optical depth to the 
observed value. 



Finally we would like to characterize the excess reddening of the microlensed stars. This can 
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be done by denning an observable parameter ^ Q b s such that 

^obs)^obs / 00 \ 
° bs ~~ JdK ohs N s (K ohs ) JdK ohs N*(K ohB ) 

JdDv s (D)K(D) _ JdDMD) K (D) 

JdDv s {D) JdDu*(D) ' 1 ' 

Because of the integration over the entire range of the reddening, £ b s is independent of the 
dispersion a between the observed reddening rank of a star and the reddening rank predicted from 
a smooth model (cf. eq. [Ll|), 



3. Results 

3.1. Model parameters 

Here we study the effect of a mixed halo model on the reddening distribution of the 
microlensed stars. Since the parameters of the LMC are very uncertain, we will explore a range of 
dust models and stellar models. Table 1 lists the parameters for the models. 

We fix the distance to the LMC -Dlmc = 50kpc, and set the tidal radius of the LMC 
l-D min — -DlmcI = l-Cmax — -Dlmc I = lOkpc. For the sech 2 stellar disk of the LMC and the dust 
layer, we fix S disk = 4OOM pc~ 2 . We fix the FWHM thickness of the dust layer w = 200 
pc, but allow the relative thickness of the two disks W/w to vary in the range 0.5 to 2, a 
reasonable range for thin stellar and dust disks. We set the FWHM of the tidal material 
WW = 0.05(D 

max — -Dmin) — lkpc, but allow the peak position to vary in the range 40 
kpc < -Dextra < 60 kpc, and the amount of stars in extra component to vary in the range 
< S oxtra < 50M Q pc~ 2 . We set this upper limit for the surface density of the extra material to 
about 10% of the surface density of the LMC, an acceptable amount for some hidden material. In 
comparison the surface density of dark matter (machos plus wimps) halo Ehalo ~ 100M Q pc~ 2 (cf. 
eq. [?]) . These parameters are comparable to those of previous models for the volume density of the 
LMC disk (Wu 1995, Weinberg 1999). 

It turns out the reddening distributions of microlensed and unlensed stars are insensitive 
to exact values of w, W, W e xtra an d Sdisk as long as we are in the regime where 
W ~ W <C Wextra ^ l^max — -D min | and 

^extra ~ ^haio ^ ^disk- These distributions are 
more sensitive to the ratios such as W/w, £ ex t ra /£iiaio! an< ^ j^ xtra ~^ LMI ?l , For example, models 

l-L-'max — i-'minl 

with the extra stars centered on the LMC disk, i.e., l^ xtra ~^ LM c 0, require much more stars 
with Scxtra/Shaio ^ 1 for come up with enough events, hence are much less efficient than models 
with extra stars 5 — lOkpc in front or behind the LMC in terms of producing microlensing. So 
the main parameters that we will vary are S ex t ra , -Dextra an d W/w. The predictions are made for 
halo-lensing models and mixed models with different amount of stars in the extra component and 
different offset distance from the LMC and the relative thickness of the stellar disk vs. the dust 
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disk. 

We also set the random error a for the observed reddening at 20%-40% level, which seems a 
reasonably amount. A larger dispersion would smooth out some tell-tale features of the reddening 
distributions. But it turns out that the mean excess £ Q b s and and the optical depth in the 
mid-plane T b s (K bs = 1) are insensitive to the observational error a. 

3.2. Reddening distributions of microlensed and unlensed stars 

The line of sight distributions of a few model quantities are shown in Fig. [I] for one of the 
models (see caption). Clearly all stellar density distributions peak at D = Dlmc = 50kpc, the 
mid-plane of the high density disk of the LMC; similarly for the thin dust layer. The LMC stars 
have a second peak due to extra stars placed behind the LMC disk. The lens density distribution 
also has a gentle falling part in the range 10 kpc < D < 40 kpc due to foreground machos, whose 
number density falls as D~ 2 . The lensed stars (dotted line with diamonds) also have a higher tail. 

For sources in front the LMC disk D < -Dlmc the microlensing optical depth is nearly 
constant, r(L> s )/r st d ~ /macho ~ 0.36. This is because the optical depth due to the ever decreasing 
density of the foreground machos is insensitive to the source distance. But after passing the 
LMC disk the optical depth takes off linearly with the source distance till T(-D s )/r s td S> 1- This 
is because any source star behind the LMC suddenly sees an intervening dense screen of lenses 
coming from the high surface density stellar disk of the LMC, so the optical depth goes up in 
proportional to the distance to the LMC disk (cf. eq. ^). This explains the higher tail of the 
lensed stars as compared to unlensed stars (cf. Fig. |l|). Note that there is only a modest increase 
of the optical depth from entering to exiting the LMC thin disk. The amount is on the order of 
(^disk/Shaio)(^/10kpc) < 8% of the total optical depth, consistent with the estimation of Sahu 
(1994) and Wu (1994) for a thin disk. 

The dust extinction is nearly a Heaviside function of the distance, climbing steeply from zero 
absorption in front of the layer to a constant value 2 x ^4lmc(-Dlmc) a few hundred pc behind 
the layer, where ^4lmc(-^lmc) is the absorption in the mid-plane. This is because as a star sinks 
deeper in the dust layer it experiences an increasing reddening before it emerges from the dust 
layer again. 

Fig. |2] shows the reddening distributions of microlensed and unlensed stars. The distribution 
for unlensed stars (upper panels) hardly depends on parameters of the extra component S extra 
and Dextra- It is only sensitive to the relative thickness of the stellar disk and the dust layer. 
At one end of the extreme, we have a thin dust layer and a thicker stellar disk (W/w = 2, solid 
lines). We see the familiar double-horn structure. This is because most stars in a thick disk 
are either in front of the dust layer and free from reddening or behind the layer and reddened 
by the maximum amount. At the opposite extreme, if the stellar disk is thinner than the dust 
layer (W/w = 0.5, dash-dot-dot-dot lines), the distribution becomes peaky in the middle because 
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most stars of the thin disk are in the mid-plane of the dust layer, hence reddened by half of the 
maximum (back-to-front) value. In between if the dust layer and the stellar disk have the same 
scale height and are evenly mixed (W/w = 1, dashed lines), we have a flat top-hat distribution. 

The distributions for the microlensed stars depend on the location and the amount of the 
extra component. For halo-lensing models, the distribution of lensed stars and unlensed stars are 
virtually indistinguishable (cf. Fig. ^a). The situation is largely unchanged if we move all lenses 
from the halo to 10 kpc in front of the LMC disk (cf. Fig. §b). The extra component shows up 
merely as a marginal excess of low reddening stars in the unlensed stars, a small effect which 
might escape detection. 

The situation is drasticly different if we put even a small amount (S extra = IOMq pc~ 2 , or 
2.5%) of stars behind the LMC disk (Fig. ||c). The distribution are strongly skewed towards high 
reddening. This is because lensing favors source stars well-behind the lenses. More specificly if 
the lenses are machos midway to the LMC, or stars some 10 kpc in front of the LMC, then the 
probability of drawing a source star at 100 pc behind the dust layer is about 0.8% or 2% higher 
than drawing one at 100 pc in front of the dust layer (cf. eq. ^). So in average, the lensed sources 
should be more reddened than the unlensed ones by merely a few percent. In contrast, if the lenses 
are at mid-plane of the LMC, then nearly all sources will be behind the dust layer, and hence the 



reddening rank K Q bs increases by one (cf. eq. || and eq. 10). 

In Fig. |2|c there are still some source stars at low reddening because we allow for machos in 
the foreground. The skewness effect becomes stronger as we decrease the macho fraction, and 
increase the extra stars behind the LMC. For X e xtra = 4OM pc -2 and -D C xt m = 60kpc, then all 
lensed sources should have an observed reddening rank K bs 

2. In this case / mac ho ~ and all 

lenses come from the LMC disk. 

Another way to quantify the above effects is to look at the observed optical depth as a 
function of the reddening rank of stars (cf. Fig. ||a). The plateau of the T b s (K bs) curves in the 
range < K Q b s < 1 are for stars imbeded in the dust layer of the LMC. They are slightly closer to 
us than the stars at the mid-plane, for which K Q b s = 1 ± 3a, but are virtually at the same distances 
-Dlmc =t w/2. The plateau also set an upper limit on the optical depth of machos (cf. eq. ^C|). 

On the other hand those stars with K b s > 2 include both stars well behind the LMC disk, 
which have a very high optical depth (cf. the curve of t(D) in Fig. |l[), and any Gaussian tail of 
the LMC disk stars due to random error of the reddening. So T b s (K bs) can have a high tail if the 
extra component behind the LMC is stronger than the Gaussian tail due to a non-zero dispersion 
a. These results depend somewhat on the reddening and the relative thickness of the stellar 
disk and the dust layer (cf. Fig. |||b). A larger measurement error smoothes the sharp transition 
from the plateau to the high tail. The tail is lower if the stellar disk is thicker than the dust 
layer, because of more overlap in the observed reddening of the stars in the LMC disk and those 
well-behind the LMC. 



Fig. |I| shows that in general a strong excess in reddening is seen in models with most of the 
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lenses in the LMC disk, and the sources behind the midplane of the LMC disk. Models with a 
high macho fraction / mac ho > 0.4 and/or with lenses being extra stars a few kpc well in front of 
the LMC disk predict only a modest amount of excess £ b s < 0.1. So these models could be ruled 
out if we measure a strong excess. In particular, if £ b s > 0.5 then the dark halo cannot have more 
than 15% in machos. On the other hand, if the excess is small, then the interpretation remains 
non-unique. These results depends only slightly on the thickness of the dust layer with the thinner 
dust layer predicting a stronger excess. Given accurate reddening, it is possible to measure or at 
least set an upper limit on the fraction of dark matter in machos. 



3.3. Empirical relations 

We have also find a few empirical relations for converting the reddening excess £ b s to the 
macho fraction / mac ho- Assuming that there are few stellar lenses in the immediate foreground of 
the LMC disk, which would be the case if 0.8 > £ bs > 0.2, then to a good approximation, the 
fraction /* of microlensing due to LMC stars is given by 

/* ~ 1.2Us. (25) 

The fraction of machos is related to £ Q b s by 

/macho ~ 0.5(1 - 1.2£ obs ) (26) 

. Here we have adopted the present reported value for the microlensing optical 
depth to the LMC, i.e., about half of that of the standard halo. But the results can be rescaled 
proportionally as more data become available. 

The result can also be recast as following, 

, _ T obs (0 < At obs < 1) , . 

/macho ~ g x 1Q _ 7 (41 ) 

(cf. eq. ^), where r o b s (0 < K bs < 1) is the optical depth of the LMC stars observed with lower 
than average reddening; these are stars sandwiched in a thin layer between the mid-plane of the 
LMC disk and the near side of the dust layer. The optical depth of the least reddened stars best 
approximates the optical depth of foreground machos, because it is least contaminated by stellar 
lenses in the LMC. 

These relations work best if we can neglect stellar lenses in the immediate foreground of the 
LMC disk. But interestingly, these results are insensitive to the assumption of the dust distribution 
and stellar distribution of the LMC, e.g., the distance of the extra background material, and their 
surface density, the relative thickness of the dust layer and the stellar disk. We also allow for a 
realistic amount of measurement error and patchiness of extinction. 



(cf. eq. E3 and 21 
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4. Summary 

In summary, we have studied effects of dust layer in the LMC on the microlensing events in 
a wide range of models of star and dust distributions in the LMC. We propose to bin LMC stars 
according to their reddening rank K Q b s (as defined in eq. ||), and study the observable microlensing 
optical depth as a function of the reddening rank K Q bs- We find that self-lensing models of the 
LMC draw preferentially sources behind the dust layer of the LMC, and hence can be distinguished 
from the macho-lensing models once the reddening by dust is measured. 

If a low excess in reddening (£ bs < 0.2) is observed then the interpretation would not be 
unique: the reddening distribution does not distinguish very well lensing by machos and lensing 
by stars a few kpc in front of the LMC. The optical depth can be explained equally well by a 
screen of macho lenses with S ex t; ra = 50M Q pc 2 at 10 kpc from the Sun, or by a screen of stellar 
lenses with a same surface density at 10 kpc in front of the LMC disk. As far as the reddening 
distributions are concerned, these two models are barely distinguishable (cf. Fig. 0ab). On the 
other hand the reddening distribution is very sensitive to any stars behind the LMC disk. If 2-3% 
of the LMC stars is put in an extra component 5 — 10 kpc behind the LMC, then we can see 
markedly different distribution in the microlensed stars (cf. Fig. ^c); this corresponds to a surface 
density S ex tra = 1OM pc -2 . The signal gets even stronger as we decrease the macho fraction and 
increase the surface density of the background stars. 

Our main finding is that among stars of different reddening rank, the optical depth of the 
least reddened stars r b s (Kobs ~ 0) is the closest approximation to the optical depth of machos (cf. 
Fig. ||). The macho fraction / mac ho also correlates tightly with the excess reddening £ Q b s of the 
microlensed sources (cf. Fig. Q). These results are summarized by the empirical relations in the 
previous section. An observable high excess in reddening (1 > £ b s > 0.2) would be a definative 
signal for many LMC disk lenses and low macho fraction. Potentially we can constrain the fraction 
of Galactic dark matter in machos this way. 

There are a number of problems to apply the method to observations, (a) Reddening of 
individual stars is difficult to measure accurately. Reddening can be determined by constructing 
reddening-free indices with photometry of three or more broad bands, or with low resolution 
spectroscopy. Typical accuracy is about 0.02 mag. in E(B — V) with these methods (Harris et al. 
1997), which is about 20% of the typical reddening in the LMC (Harris et al. 1997). (b) Stars in 
the LMC are likely crowded from the ground. Blended images of stars lead to unphysical colors 
and spurious reddening, (c) Dust distribution is very clumpy. Extinction towards OB stars in the 
LMC disk can easily vary at a factor of two level among different patches of the sky separated by 
l' (Harris et al. 1997). (d) The effect of extinction by the dust layer of the Galaxy near the Sun 
should be included. Oestriecher et al. (1996) show that the Galactic foreground extinction is not 
entirely smooth. There are dark patches on 30' scales. 

The dumpiness of the dust, together with the fairly large error of the reddening vector 
measurable from broad-band photometry, can lead to a large scatter in the relation between 
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reddening and line of sight depth. Zhao (1999c) argue that these problems, particularly patchiness 
can be overcome if we restrict to the smallest patches of sky around each microlensing source 
where the variation is likely at only 20% level. Spectroscopic observation with the Hubble Space 
Telescope in the ultra-violet can in principle resolve many faint stars in the small patches around 
each microlensed stars, and measure the reddening accurately. By obtaining the reddening of 
stars in the present 20 — 30 microlensing lines of sight we should be able to measure the excess 
reddening confidently. We should also be able to bin the events in the reddening rank, and study 
the optical depth as a function of reddening. A null result of any variation of the optical depth 
would be in favor of a significant baryonic dark component of the Galaxy. 
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Table 1. Fixed parameters for the stellar, dust and halo models 



^disk 


Wextra 


-^min 




w 


(Tabs) 


Tstd 


^halo 


4OOM pc" 2 


lkpc 


40kpc 


60kpc 


200pc 


2.9 x 10~ 7 


4.7 x 10~ 7 


2OOM pc" 2 
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Table 2. Variable parameters for the stellar, dust and halo models 



^extra -^extra W/ W O /macho 



O-5OM pc" 2 40-60 kpc 0.5-2 20%-40% 0-0.5 
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Fig. 1. — Upper panel: the density (in logarithmic scale with arbitrary zero point) of the lenses, 
the unlensed LMC stars and the density of the dust layer along the line of sight to the LMC. 
Middle panel: a zoom-in of the upper panel for the unlensed LMC stars and the lensed source 
stars. Lower panel: the run of the rescaled absorption of the dust (dashed lines with error bars) 
and the rescaled optical depth of microlensing (dashed line with squares) as functions of line of 
sight distance. The small error bars indicate a 20% dispersion due to patchiness and measurement 
error. All calculations are done for W = 2w = 400pc, Astra = 53kpc, S ex tra = 2OM pc -2 and 

/macho = 0.38. 
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Fig. 2. — the frequency of finding stars with a certain rescaled reddening for the microlensed 
sources (the lower part of each panel) and the unlensed stars (the upper part of each panel). The 
legends show model parameters. Models are calculated for a thick disk and a thin dust layer with 
W = 2w = 400pc (solid lines), an equally thin disk and dust layer with W = w = 200pc (dashed 
lines), and a thin disk and a thick dust layer with W = w/2 = lOOpc (dash-dot-dot-dot lines). 
The distributions have been smoothed by a Gaussian with a dispersion a = 20% to allow for the 
measurement error and patchiness. 
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Fig. 3. — the optical depth, rescaled by that of the standard model r st d, as a function of the 
reddening for several models with S cxtra = 50M Q pc~ 2 . Models are calculated for a = 20%, 
W = w = 200pc, and a range of -D ex tra = 41 — 57 kpc (left panel), and for D extia = 57 kpc, and a 
range of a = 20% — 40% and W/w = 0.5 — 2 (right panel). In the absence of stellar lenses, a halo 
model would predict the (horizontal dotted) line r bs/7"std = 0.5. 
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Fig. 4. — the excess reddening vs. the fraction of machos. We vary surface density S ex tra = 
(10 — 5O)M pc" 2 of the extra stars and its peak position D extra = 40 — 60 kpc. Each model allows 
a range of thickness for the dust layer w stellar disk and W. We trade between the amount of 
extra stars S ex tra and the amount of machos / mac ho such that models all have a fixed optical depth 
(fobs) Astd = 0-5- The predictions are shown as three hatched bands. 



